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Abstract. We discuss the efficiency of stellar gravity torques as a mechanism to account for the feeding of the central engines 
of four low luminosity Active Galactic Nuclei (AGN): NGC4321 (HII nucleus/LINER), NGC4826 (HII nucleus/LINER), 
NGC4579 (LINER 1.9/Seyfert 1.9) and NGC6951 (Seyfert 2). These galaxies have been observed as part of the NUclei 
of GA/axiei-(NUGA) CO project, aimed at the study of AGN fueling mechanisms. Our calculations allow us to derive the 
characteristic time-scales for gas flows and discuss whether torques from the stellar potentials are efficient enough to drain the 
gas angular momentum in the inner 1 kpc of these galaxies. The stellar potentials are derived using high-resolution near infrared 
(NIR) images and the averaged effective torques on the gas are estimated using the high-resolution (~0.5"-2") CO maps of the 
galaxies. Results indicate paradoxically that feeding should be thwarted close to the AGNs: in the four cases analyzed, gravity 
torques are mostly positive inside r~200 pc, resulting in no inflow on these scales. As a possible solution for the paradox, we 
speculate that the agent responsible for driving inflow to still smaller radii is transient and thus presently absent in the stellar 
potential. Alternatively, the gravity torque barrier associated with the Inner Lindblad Resonance of the bars in these galaxies 
could be overcome by other mechanisms that become competitive in due time against gravity torques. In particular, we estimate 
on a case-by-case basis the efficiency of viscosity versus gravity torques to drive AGN fueling. We find that viscosity can 
counteract moderate-to-low gravity torques on the gas if it acts on a nuclear ring of high gas surface density contrast and ~a 
few 100 pc size. 

We propose an evolutionary scenario in which gravity torques and viscosity act in concert to produce recurrent episodes of 
activity during the typical lifetime of any galaxy. In this scenario the recurrence of activity in galaxies is indirectly related to 
that of the bar instabilities although the active phases are not necessarily coincident with the maximum strength of a single bar 
episode. The general implications of these results for the current understanding of fueling of low-luminosity AGN are discussed. 

Key words. Galaxies:individual:(NGC4321, NGC4579, NGC4826, NGC6951) - Galaxies:ISM - Galaxies:kinematics and 
dynamics - Galaxies:nuclei - Galaxies:Seyfert - Radio lines: galaxies 



1. Introduction 

The phenomenon of nuclear activity is understood to be a re- 
sult of the feeding of supermassive black holes (SMBHs) in 
galactic nuclei. Observational evidence accumulated over the 
last decade indicates that SMBH exist in most galactic bulges 
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(e.g., Kormendy & Richstone I1995I Magorrian et al. [1998; 
Ferrarese & Merritt 2000, Gebhardt et al. I2000> . Among all 
these massive black holes, very few are highly active. AGN 
are found in 10% of the local galaxies (Ho et al. 119971 ); how- 
ever, this percentage is increased up to ~44% if LINERs are 
taken into account. One of the present challenges is to under- 
stand how AGN can be fed during their lifetime. In the feed- 
ing problem the gas supply must come from the whole disk 
of the host galaxy at large distances compared to the radius 
of gravitational influence of the central engine. Therefore it is 
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expected that a hierarchy of mechanisms combine to drive vir- 
tually all the gas from the large ~kpc scales down to the inner 
~pc scales. The different spatial scales involved suggest that 
the various mechanisms at work have very different time scales 
(Shlosman et al. rE%9IIT9901 Combes l^UOlllSUOll Jogee EOOfl . 
Recent observational and theoretical evidence indicates that the 
AGN lifetimes may be as short as ~a few 10 7 -10 8 yr (Ho et 
al. l2003l Martini l2004l Wada l2004l Merloni l2004l ). Moreover, 
Wada (2004) finds evidence that mass accretion may not be 
constant even during the nominal duty cycle of 10 8 yr, but com- 
posed of several shorter episodes with a duration of 10 4_5 yr. 
This time-scale conspiracy could explain the lack of success 
of observers in finding any correlation between the presence 
of ~kpc scale non-axisymmetric perturbations (e.g., large- 
scale bars and interactions) and the onset of activity in galax- 
ies, except for very high luminosity objects(QSOs) (Moles 
et al. [T9931 Mulchaey & Regan \TWH Knapen et al. 2000 
Krongold et al. 120011 Schmitt 2001). On these spatial scales, 
Hunt & Malkan ( 1999) have found a significant correlation be- 
tween the detection rate of outer rings and the onset of activity. 

The search for a universal feeding mechanism has been 
pursued by looking for morphological features in the cen- 
tral kpc of nearby AGN with high spatial resolution (~a few 
100 pc), though with limited success. Nuclear stellar bars 
seem to be as common in AGN as in non-AGN (Regan & 
Mulchaev l 19991 Laine et al.[2002i. Furthermore, a similar sce- 
nario holds for nuclear spirals: while Martini & Pogge ( 1999) 
and Pogge & Martini ( 2002 ) initially found a high frequency of 
dusty spirals in the HST enhanced color images of Seyfert nu- 
clei of their sample, more recently, Martini et al. (J3003 ) have 
shown on a firmer statistical basis that these nuclear features 
are not preferentially found in AGN. Similarly there is only 
weak statistical evidence that nuclear rings are more frequently 
found in Seyferts (Knapen l2"0"05V Hunt & Malkan ( 2004 ) found 
evidence for an excess in type 2 Seyferts of kpc-scale twisted 
isophotes, though the cause of such twists is unclear. 

On the modeling front, significant progress has been made 
in recent years on the study of the feeding efficiency of dif- 
ferent types of gravitational instabilities: nested bars (e.g., 
Shlosman et al. 119891 Friedli & Martinet 1993 Maciejewski 
& Sparke 120001 Englmaier et al. 2004), gas spiral waves 
(e.g., Englmaier & Shlosman 2000 Maciejewski et al. 120021 
Maciejewski 2004a 2004b), m - 1 perturbations (e.g., Shu 
et al. 119901 Junqueira & Combes 119961 Garcfa-Burillo et 
al. |2000) and nuclear warps (e.g., Schinnerer et al. 2000). 
However, the lack of high quality multi-wavelength observa- 
tional constraints on the different models has thus far made the 
choice of a single optimal scenario rather difficult. 

The study of interstellar gas in AGN is essential to under- 
stand the phenomenon of nuclear activity in galaxies and its 
possible link to circumnuclear star formation. As most of the 
neutral gas in galactic nuclei is in the molecular phase, CO 
lines are best suited to undertake high-resolution mapping of 
AGN hosts, with interferometer resolution of <100 pc, i.e., 
the scales on which secondary modes embedded in kpc-scale 
perturbations are expected to take over. CO lines better trace 
the total gas column densities than dust extinction probes ob- 
tained from HST NIR/optical color images. Most importantly, 



CO maps provide the gas kinematics (velocity fields and ve- 
locity dispersions). This information is essential to charac- 
terize gravitational instabilities and to constrain the models. 
The NUclei of GAlaxies-NUGA-project, fully described by 
Garcfa-Burillo et al. (2003a 2003b), is the first high-resolution 
(~0.5"-l") CO survey of 12 low luminosity AGN (LLAGN) 
including the full sequence of activity types (Seyferts, LINERs 
and transition objects from HII to LINER). In the case of 
LLAGN, the required mass accretion rates derived from the 
typical bolometric luminosities of these objects range from 
10" 2 to 10 _5 M G yr _1 (from Seyferts to LINERs; e.g., see compi- 
lation by Jogee 2004 ). Observations, carried out with the IRAM 
Plateau de Bure Interferometer (PdBI), have been completed 
early 2004. NUGA surpasses in both spatial resolution and sen- 
sitivity ongoing surveys of nearby AGN conducted at OVRO 
(MAIN: Jogee et al. l20T)Tl and at NRO (Kohno et al. l200H . 

In this paper we focus on the study of gravitational torques 
in a subset of NUGA galaxies, which span the range of 
the different activity classes within our sample: NGC4321 
(transition object: HII/LINER), NGC4826 (transition ob- 
ject: HII/LINER), NGC4579 (LINER 1.9/Seyfert 1.9) and 
NGC 695 1 (Seyfert 2). Information on the stellar potentials, ob- 
tained through available HST and ground-based optical/NIR 
images of the sample, is used to determine the gravitational 
torques exerted by the derived stellar potentials on the gaseous 
disk. The efficiency with which gravitational torques drain the 
angular momentum of the gas depends first on the strength of 
the non-axisymmetric perturbations of the potential (m>0) but, 
also, on the existence of significant phase shifts between the gas 
and the stellar distributions. The estimate of these phase shifts 
necessarily requires the availability of images of comparably 
high spatial resolution (<0.5" in our case) showing the distribu- 
tion of the stars and the gas. In this paper we purposely neglect 
the role of gas self-gravity as a source of non-axisymmetry in 
the gravitational potential. Nuclear galaxy disks with a high gas 
surface density and a mostly axisymmetric stellar potential can 
be prone to develop this kind of gas self-gravitating perturba- 
tion. We purposely defer the study of pure gas instabilities and 
their ability to drive gas inflow to a forthcoming publication. 

We describe in Sect. EJ the observations used, includ- 
ing high-resolution CO maps and NIR images of NGC 4321, 
NGC 4826, NGC 4579 and NGC 6951. Sect.EJinterprets them 
in terms of AGN feeding. Sect. |4] computes from NIR images 
the gravitational potentials and forces, and deduces from the 
CO maps the effective torques applied to the gas. From these 
torques, it is possible to derive time-scales for gas flows and 
discuss whether gravity torques alone are efficient enough to 
feed the AGNs. The general implications of these results for the 
current understanding of AGN feeding are presented in Sect. [5] 

2. Observations 

2.1. CO NUGA observations 

Observations of the circumnuclear disks of NGC 4826, 
NGC 4579 and NGC 6951 were carried out as part of the 
NUGA survey with the PdBI between December 2000 and 
March 2003. We used the ABCD set of configurations of the 
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Table 1. We list the position angle (PA), inclination (z), distance (D), spatial scale, coordinates of the dynamical center derived 
from CO and/or from the locus of the radio continuum source [(RA2000, Dec2ooo); accurate to ~0.5" on average] and systemic 
velocity(v^/) derived from CO for NGC 432 1 , NGC 4826, NGC 4579 and NGC 695 1 used in this paper. The spatial resolutions 
of the CO observations are listed in columns 9-10. 
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Fig. 1. a) The 12 CO(1-0) integrated intensity map obtained with the PdBI (contour levels from 1.5 to 5 in steps of 
0.5 Jy km s -1 beam -1 and from 5 to 10 in steps of 1 Jy km s -1 beam -1 ) observed in the nucleus of NGC 4321. The filled 
ellipse at the top right corner represents the CO beam size. The central r~l kpc molecular disk of NGC 4321 consists of two 
nuclear spiral arms connected (e.g., B component) to a marginally resolved r~150 pc disk of ~1O 8 M of molecular gas cen- 
tered on the AGN locus (highlighted by the star marker). {Aa, A<S)-offsets are with respect to the location of the AGN [(RA2000, 
Dec2ooo)=(12 /, 22 m 54.91' s , 15 d 49'"19.9 s )]. This position coincides within the errors with a secondary radio continuum maximum 
measured at 6 cm by Weiler et al. ( 198 1 ) and with a peak in the (J,H,K)-2MASS image of the galaxy (Jarrett et al. 2003 ). Fig. lb 
illustrates the kinematics of molecular gas along the strip at PA=127° (thick line in Fig. la). Levels go from 0.015 to 0.14 in 
steps of 0.0125 Jy beam -1 . Velocities are relative to the systemic velocity (v^/=1573 km s -1 ; determined from CO kinematics 



by Garcfa-Burillo et al. 1998 ) and Ax offsets are relative to the AGN. The major axis orientation is shown by the dashed line. Gas 
kinematics over the emission bridge B depart from circular rotation. 



array (Guilloteau et al. 1992). This assures high spatial resolu- 
tion (<1" at the highest frequency) but also an optimum sensi- 
tivity to all spatial frequencies in the maps. A previous set of 
observations of NGC 4826, using data taken with the BCD con- 
figurations, has been discussed in Garcfa-Burillo et al. (2003b) 
(hereafter called paper I). We have observed simultaneously 
the J=l-0 and J=2-l lines of 12 CO in single fields. Table □ 
lists the spatial resolution of these CO observations as well 
as other source related parameters. The primary beam size is 
42" (21") in all the 1-0 (2-1) line observations. As the bulk of 
the relevant nuclear disk emission arises well inside the central 



15" in the four galaxies discussed in this paper, observations 
have not been corrected for primary beam attenuation. During 
the observations the spectral correlator was split in two halves 
centered at the transition rest frequencies corrected for the as- 
sumed recession velocities. The correlator configuration covers 
a bandwidth of 580 MHz for each line, using four 160 MHz- 
wide units; this is equivalent to 1510km s -1 (755km s -1 ) at 
115 GHz (230 GHz). Visibilities were obtained using on-source 
integration times of 20 minutes framed by short (~ 2min) 
phase and amplitude calibrations on nearby quasars. The abso- 
lute flux scale in our maps was derived to a 10% accuracy based 
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Fig. 2. a) The 12 CO(2-l) integrated intensity map obtained with the PdBI (contour levels from 2.8 to 10 in steps of 
0.9Jy km s _1 beam 1 ) is overlaid on the B-I color image from HST (grey scale) observed in the nucleus of NGC4826. The 
filled ellipse at the top right corner represents the CO beam size. (Aa, A<5)-offsets are with respect to the location of the AGN 
(marked by the star): (RA 2 ooo, Dec 2 ooo)=(12 /l 56 m 43.63 s , 21 c, 40 m 59.1 s ). The AGN locus is identified by a blue point-like source 
in the B-I map; it also coincides with a non-thermal radio continuum peak measured at 6 cm by Turner & Ho ( 1994 1. The 80 pc 
radius circumnuclear disk (CND) of NGC4826, with ~3xlO 7 M of molecular gas, shows a lopsided ringed disk morphology; 
the disk is off-center with respect to the AGN. b) The kinematics of the gas, here displayed along the minor-axis p-v plot (thick 
line in Fig. 2a), are suggestive of strong streaming motions at the crossing of the ring edges (N, S). Levels go from 0.025 to 0. 19 
in steps of 0.015 Jy beam -1 . Velocities are relative to the systemic velocity (v^ s R =413 km s ; determined from CO kinematics 
by Garcfa-Burillo et al. 2003b I and Ay offsets are relative to the AGN. 



on the observations of primary calibrators whose fluxes were 
determined from a combined set of measurements obtained at 
the 30m telescope and the PdBI array. Image reconstruction 
was done using standard IRAM/GAG software (Guilloteau & 
Lucas 2000 1. In this work we use prior observations of the 1— 
line of 12 CO of NGC4321 made using the BCD configura- 
tions of the PdBI and previously published by Garcfa-Burillo 
et al. dl998> (see this paper and Table^for details). 

Depending on the resolution/sensitivity requirements, we 
use either naturally or uniformly weighted line maps, as in- 
dicated throughout the paper. Uniform weighting in the 2-1 
line enables us to achieve subarcsecond spatial resolution in 
the maps of NGC 4826, NGC 4579 and NGC 695 1 . By default, 
all velocities are referred to the systemic velocities (v jyj , listed 
in Table as determined from this work and from Garcfa- 
Burillo et al. ( 1998l. Similarly (Aa, A6) offsets are relative to 
the AGN loci derived from our own estimates (Table [Q. 

Molecular gas masses are derived from the CO(1-0) in- 
tegrated intensities assuming a CO-to-H2 conversion factor 
X=N(H 2 )/Ico(i-0)=2.2xl0 20 cm- 2 KT 1 kirr 1 s (Solomon & 



Barrett 1991 1. When required, molecular gas column densities 
are inferred from the CO(2-l) integrated intensity maps. In this 
case CO(2-l) intensities are first corrected by the 2-1/1-0 ra- 
tio measured within the equivalent 1-0 beam at each position, 
and subsequently we convert them into N(H 2 ) assuming the X 
factor referred to above. 

We have estimated the percentage of CO(1-0) flux recov- 
ered in the PdBI maps by comparing the single-point fluxes 
detected by the 30m telescope (Garcfa-Burillo & Krips, private 
communication) towards the nuclei of the four galaxies dis- 
cussed in this paper with the fluxes recovered in the PdBI maps, 
corrected by primary beam attenuation and convolved to the 
30m resolution at this frequency (21"). The fraction of the flux 
recovered inside the 21" field-of-view ranges from 60%-65% 
in NGC 6951 and NGC 4579 to 75% in NGC 4321 and 90% 
in NGC 4826. The corresponding values for the 2-1 line of CO 
are comparable. The missing zero spacing flux in these maps is 
expected to be found in low-level emission arising in the shape 
of smooth extended components. As has been shown by Heifer 
et al. ( 2003 1, who estimated the percentage of flux recovery in 
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Fig. 3. a) The 12 CO(2-l) integrated intensity map obtained with the PdBI (contour levels from 0.4, 0.7, 1.3, 2.0 to 1 1 in steps of 
0.9Jy km s _1 beam 1 ) is overlaid on the V-I color image from HST (grey scale) observed in the nucleus of NGC4579. The filled 
ellipse at the top right corner represents the CO beam size. (Aa, A(5)-offsets are with respect to the location of the AGN (marked 
by the star): (RA 20 oo, Dec 2 ooo)=(12 /! 37 m 43.52 ? , ll' , 49 m 05.5 s ). The AGN locus is identified by a point-like continuum source 
detected at 3 mm and 1 mm by Garcfa-Burillo et al. 2005 (in prep). Two spiral arcs concentrate the bulk of the ~3.2xl0 8 M Q 
molecular gas mass in the central 1 kpc of the galaxy. A close-up view of the inner 200 pc region (shown in c), shows a central 
ringed disk (highlighted by the dashed ellipse) with the AGN lying on its southwestern edge. A gas clump of ~10 6 M o (denoted as 
E) delimits the disk to the East, b) The kinematics along the major axis (thick line in Fig. 3a) reveal highly non-circular motions 
related to the E clump. Levels go from 0.006 to 0.09 in steps of 0.006 Jy beam -1 . Velocities are relative to the systemic velocity 
(Vj V ^ R =1469 km s _1 ; determined from CO kinematics by Garcfa-Burillo et al. 2005, in prep.) and Ax offsets are relative to the 
AGN. d) Same as b) but here along the declination axis (dashed thick line in Fig. 3a) with levels going from 0.006 to 0.024 in 
steps of 0.006 Jy beam -1 . Highly non-circular motions are related to the N clump. 



BIMA SONG galaxies using the 12m NRAO telescope, this 
percentage is always very high in the central regions of galax- 
ies (i.e., the domain of NUGA maps). The reason is that the ve- 
locity gradient is largest at the nuclear regions. CO emission is 
thus confined to much smaller areas in individual channel maps 
compared to the outer disk regions where the percentage of flux 
filtered out can be higher. Considering that the percentage of 
flux actually present in the PdBI maps over the total single- 



dish estimate is moderate-to-large in the galaxies studied here 
(see above), we do not expect that the morphology of the maps 
will significantly change by the addition of a plateau-like com- 
ponent. Moreover, the gravity torque calculation developed in 
this paper is based on azimuthal averages made to infer time- 
scales for the gas flows in these galaxy nuclei. The approach 
followed makes our results virtually insensitive to the presence 
of a weak extended component: gravity torques on this type of 
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Fig. 4. a) The 12 CO(2-l) integrated intensity map obtained with the PdBI (contour levels from 0.15, 0.25, 0.40, 0.70 to 6.2 in 
steps of 0.5Jy km s beam 1 ) is overlaid on the J-H color image from HST (grey scale) observed in the nucleus of NGC 695 1 . 
The filled ellipse at the top right corner represents the CO beam size. (Aa, A<5)-offsets are with respect to the location of the AGN 
(marked by the star): (RA 20 oo, Dec 2 ooo)=(20 /l 37 m 14.12 s , 66 rf 06'"20.0 5 ). The position of the AGN is given by the point-like radio 



continuum source measured at 6 cm and 20 cm by Ho & Ulvestad (2001 1. The molecular gas distribution in the central 1 kpc 
shows two highly contrasted nuclear spiral arms containing 3xlO 8 M which are presently feeding a circumnuclear starburst. A 
compact molecular complex (denoted as G) of ~a few 1O 6 M is detected on the AGN. Furthermore, we have tentatively detected 
a northern molecular gas component linking G with the N spiral arm (i.e., the spiral running North from West) which could be 
related to the filamentary dusty spiral seen in the J-H color HST image (see also Fig. 4b). b) The kinematics along the major 
axis (thick line in Fig. 4a) are compatible with circular motions for the gas near the AGN. Levels go from 0.008, 0.0012, 0.016 
to 0.076 in steps of 0.006 Jy beam -1 . Velocities are relative to the systemic velocity (v" R =1441 km s -1 ; determined from CO 
kinematics by Schinnerer et al. 2005, in prep.) and Ax offsets are relative to the AGN. 



source distribution will be zero. Therefore we do not expect the 
derived torque budget to be significantly biased. 

2.2. Near-infrared and optical observations 

We acquired from the HST archive 1 broadband images 
of NGC 4826, NGC 4579 and NGC 6951, including three 



1 Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Science Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 5-26555. 



NICMOS images (F110W and F160W for NGC 6951; F160W 
for NGC 4826) and four WFPC2 images (F450W and F814W 
for NGC 4826; F555W and F814W for NGC 4579). The op- 
tical images were combined using (crreject) to eliminate cos- 
mic rays, and calibrated according to Holtzman et al. {T995 1. 
The NICMOS images were re-reduced with the STSDAS task 
calnica using the best reference files, and the images were cal- 
ibrated in the standard way. The "pedestal" effect (see Boker 
et al. 119991 was removed with the van der Marel algorithm 2 . 
Sky values were assumed to be zero since the galaxy filled 



2 http://www. stsci. eduf marel/software/pedestal. html 
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the WFPC2/NICMOS frames, an assumption which makes 
an error of ~0.1 mag at most, in the corner of the im- 
ages. A J-H color image of NGC 6951 was constructed from 
Fl 10W-F160W according to the transformations by Origlia & 
Leitherer (2000|. For NGC 4321, we have adopted the (ground- 
based) K-band image presented by Knapen et al. d!995> . 



3. Observational evidence of ongoing feeding 

3.1. NGC 4321 

The inner r~1.5 kpc of this galaxy was mapped by Garcfa- 
Burillo et al. d!998i with the PdBI at moderate (~2") spatial 
resolution in the 1-0 line emission of 12 CO. NGC 4321 has 
been classified as a transition object (i.e., HII/LINER) by Ho 
et al. ( 119971 1. As can be seen in Fig. molecular gas in the 
nucleus of NGC 4321 is concentrated in a two spiral arm struc- 
ture that starts at r~550 pc, near the end points of a prominent 
nuclear bar (detected in the K band by Knapen et al. 1995J see 
also Fig.|5jl and extends out to r~1.2 kpc. There is also a cen- 
tral CO source coinciding with the AGN which is marginally 
resolved (r~150 pc) by the ~2" beam. This central source con- 
tains a molecular gas mass of ~1O 8 M . The CO spiral arms 
mostly lie at the trailing edges of the nuclear bar. This particu- 
lar geometry determines the feeding budget for the gas in this 
region (see Sect. 14.2.11 . The spiral arms and the central source 
are connected by a molecular gas bridge which is spatially re- 
solved NW of the nucleus, i.e., at the leading edge of the nu- 
clear bar (component B in Fig. [2?)- Garcia- Burilloet al. (1998) 
interpreted these CO observations and their relation with other 
gaseous and stellar tracers using numerical simulations of the 
cloud hydrodynamics. They found that the best fit for the gas 
flow corresponds to the nuclear bar being decoupled from the 
large-scale bar. The nuclear bar is fast, with a pattern speed 
Sl p >150 km s _1 kpc _1 . This pushes corotation of the nuclear 
bar inward, likely inside the outer edges of the CO spiral arms. 

The kinematics of molecular gas are characterized by 
streaming motions detected in the CO spiral arms (see dis- 
cussion in Garcfa-Burillo et al. 119981 . Closer to the AGN 
(r<500 pc), we find that molecular gas also displays significant 
departures from circular rotation at the location of B. Fig. ^> 
shows the position-velocity (p-v) plot along a strip at PA=127°, 
purposely oriented to illustrate the CO kinematics at B. The 
discontinuity in radial velocities at Ax~-4" indicates that the 
molecular gas flow is decelerated, exactly as expected if gas 
flows along the leading edges of the bar. Details on the kine- 
matics of molecular gas in the central component are hidden 
due to the insufficient spatial resolution of the CO maps. We 
notice however that molecular gas emission is asymmetric with 
respect to v sys : CO emission is preferentially blue-shifted. 

As is the case for the other transition object analyzed in this 
paper (NGC 4826), we have found a large molecular gas con- 
centration (~1O 8 M ) near the AGN in NGC 4321 (r<150pc). 
The analysis of gas kinematics provides evidence for gas fuel- 
ing at present on intermediate scales: at r~200-300pc from the 
AGN. However, the spatial resolution of the CO maps does not 
allow us to probe closer than 100 pc from the AGN. 



3.2. NGC 4826 

The first NUGA maps of the transition object NGC 4826 
were published in paper I. The CO images showed already a 
large concentration of molecular gas (~16OOM0/pc 2 ) in the 
~160 pc-diameter circumnuclear disk (CND) of this galaxy. 
The distribution of molecular gas in the inner CND is signif- 
icantly lopsided with respect to the position of the AGN; this 
suggests that m=l instabilities may be at work at radial dis- 
tances of r~50-60 pc from the central engine. With the newest 
0.5" (10 pc) resolution 12 CO(2-l) observations the distribu- 
tion of molecular gas in the CND is fully resolved; the CND 
appears in the new maps as an off-center ringed disk (Fig.[2jj). 
The dynamical center of the galaxy determined from these ob- 
servations coincides within the errors with a blue point source 
identified in the B-I HST color map of Fig. [2?. This confirms 
our earlier findings that the putative super massive black hole 
lies on the southeastern inner side of the off-center ringed disk. 

The gas kinematics in the CND are characterized by the 
presence of streaming motions. A first analysis of the 2D kine- 
matics performed on the lowest resolution data of paper I in- 
dicated that the instabilities identified in the CND (and those 
of the inner m=l spiral) may not favor AGN feeding. The in- 
formation contained in the newest images confirms this result. 
Fig. [2? shows the departures from circular motion of gas ve- 
locities, identified at the northern and southern crossings of the 
CND ring along the minor axis. Gas velocities become sys- 
tematically redder (bluer) when we approach the nucleus from 
the southern (northern) side of the ring along the minor axis. 
Deprojected onto the galaxy plane (North is the near side), this 
pattern indicates that the radial velocity component changes 
from an inflow signature (outside the ring) into an outflow sig- 
nature (inside the ring). As it is fully discussed in paper I, this 
measured change of sign across the minor axis is compatible 
with the pattern expected for a trailing wave outside corotation 
(fast wave), i.e., the type of perturbation that would not help 
to drain the gas angular momentum. The driving agent of these 
'mainly' gaseous m = 1 instabilities may not be related to the 
stellar potential which is essentially featureless and mostly ax- 
isymmetric (see Sect. l4.2.2l and discussion of paper I). 

Based on the CO maps of NGC 4826 we find little evidence 
of ongoing AGN feeding at scales r<150pc from the AGN. 
While the molecular gas reservoir of NGC 4826 is abundant 
(>3xlO 7 M ) close to its central engine (r<150 pc), the analy- 
sis of gas kinematics provides no evidence that ongoing AGN 
feeding is at work in this transition object. 

3.3. NGC 4579 

The I2 CO(2-l) emission in the S1.9/L1.9 galaxy NGC 4579 
has been mapped at ~0.5" resolution (see Fig. 03 adapted from 
Garcfa-Burillo et al. 2005 in prep.). The molecular gas distri- 
bution in the central ~1 kpc of NGC 4579 suggests that the 
gas flow responds to the 9 kpc-diameter stellar bar identified 
in all the NIR images of this galaxy (e.g. Jarret et al. 2003 1. 
The 12 CO(2-l) map of Fig.fJJz reveals a mass of 3.2xl0 8 M o 
of molecular gas piled up in two highly contrasted spiral arcs. 
The CO lanes lie at the leading edges of the stellar bar, which is 
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oriented along PA=58°. The northern spiral is more continuous 
and better delineated than its southern counterpart; it is also 
very well correlated with the red lane seen North in the V-I 
color HST image of the galaxy, shown in Fig.[3J?. There is lit- 
tle molecular gas at r<100 pc distance from the central engine 
of NGC4579: the closest gas complex of ~1O 6 M lies East of 
the nucleus at r~150pc (complex E in Fig. |3ji, b). There is no 
molecular gas emission coincident with the AGN itself to a 3<x 
detection limit of ~a few 1O 5 M . The V-I color HST image of 
the galaxy helps to identify the eastern molecular complex as 
part of a structured disk of 150pc-diameter and ring-like shape 
(Fig. 02, c). The position of the AGN is well defined by its radio 
continuum emission detected at both 1mm and 3mm coming 
from a point source that lies close to the southwestern edge of 
the central disk. This indicates that the m — 2 point-symmetry 
of the gas flow driven by the bar of NGC 4579 breaks up at 
r<200pc and lets lopsidedness take over. An independent con- 
firmation of this picture comes from the new HST image of 
the nuclear region of NGC 4579 obtained with the ACS camera 
at 3300 A; this image resolves the central disk into a winding 
m = 1 spiral instability that mimics a ring (Contini 2004). 

The kinematics of molecular gas in the central 1 kpc of 
NGC 4579 are characterized by the presence of highly non- 
circular motions detected over the spiral arms and most no- 
tably over the central disk. Fig. [3J> shows the p-v plot along the 
kinematic major axis of NGC 4579 (PA=95°). The radial ve- 
locities of the gas in the central gas disk depart by >100kms~' 
from the expected pattern of circular rotation: emission of the 
Eastern gas complex (at Ax~1.5") appears at highly forbidden 
negative velocities (i.e., v<v svs ). Assuming that the gas flows 
inside the galaxy plane, the reported velocity deviations mea- 
sured at E would imply that gas is apparently counter-rotating 
at a speed of v~150kms _1 . As it is discussed in Garcfa-Burillo 
et al. 2005, this could be qualitatively explained by very eccen- 
tric m = 1 orbits. Alternatively, this velocity pattern could be 
accounted for assuming that gas is flowing out of the galaxy 
plane, possibly entrained by an expanding shell. The expand- 
ing shell scenario is supported by the observed kinematics of 
the gas close to the N side of the central disk. Fig. |3j/ shows 
the p-v plot taken along the declination axis, i.e., very close to 
the orientation of the minor axis. Molecular gas kinematics at 
the N complex (where v-v sys <-75 kms -1 ) can be interpreted 
either as gas flowing outward inside the plane along the minor 
axis (North is the near side), or as a signature of out of the plane 
motions (similarly to the case of the E complex). In either case 
this implies that AGN fueling is presently thwarted on these 
scales. 

In summary, most of the molecular gas content of the cen- 
tral 1 kpc of NGC 4579 is trapped in a two arm spiral struc- 
ture that can be traced from r~ 1 kpc down to r~200 pc. Some 
molecular gas (1O 6 M ) is detected at r<150 pc from the cen- 
tral engine of NGC 4579, but not on the position of the central 
engine itself (<a few 1O 5 M ). The first-order interpretation of 
the complex gas kinematics at r<150 pc provides no evidence 
of ongoing inflow from these scales down to the AGN, but on 
the contrary, it indicates outflow motions. 



3.4. NGC 6951 

NGC 6951 is a prototypical Seyfert 2 galaxy for which sub- 
arcsecond resolution 12 CO(2-l) maps have been completed 
within the NUGA project (Fig.|4j? adapted from Schinnerer et 
al. 2005, in prep.). The molecular gas distribution in the cen- 
tral 1 kpc consists of two nuclear spiral arms that can be traced 
for over 180°; the winding spiral arms end up as a highly con- 
trasted pseudo-ring at r~350 pc. The spiral arms can be iden- 
tified by their red color in the J-H HST image of the galaxy, 
shown in Fig. @J?. This m = 2 gas instability contains a sig- 
nificant gas reservoir of 3xlO 8 M which is presently feeding a 
nuclear starburst, also identified by its intense radio continuum 
and Ha emissions (Ho & Ulvestad 2001, Rozas et al. 2002). 
The geometry of the molecular gas ridges likely reflects the 
crowding of molecular clouds along the X2 family of orbits 
of the prominent stellar bar, detected in all NIR images of 
NGC 6951 ( Marquez & Moles [Tff^l Friedli et al. [19%! Perez 
et al. 120001 . As shown in Fig. |4jz, only a small amount of 
molecular gas has succeeded in making its way down to the 
AGN: most of the molecular gas mass is trapped in the nuclear 
spiral arms, quite similar to the case of NGC 4579 (though in 
this galaxy molecular gas likely populates both Xi and X2 or- 
bits). A compact unresolved molecular complex (denoted as G 
in Fig. 13?, b) of ~a few 1O 6 M is detected at the position of the 
central engine. This component could correspond to a molecu- 
lar torus (of ~40-50 pc size). As seen in Fig.^j, low level CO 
emission has been tentatively detected inside the ring, bridging 
the 'apparent' gap between the N spiral arm (running North 
from West) and the central source. The bridge is better identi- 
fied in the p-v plot of Fig.|4j?. The molecular gas mass of this 
emission bridge is more accurately estimated using the natural 
weighted map: the gas mass amounts to ~1O 7 M (Schinnerer et 
al. 2005). This could be the northern molecular counterpart of 
the filamentary spiral structure identified in the J-H HST map 

(Fig-Eb). 

The kinematics of molecular gas in the nuclear spiral 
arms reveal streaming motions (Schinnerer et al. 2005) also 
identified in the previous lower resolution CO and HCN 
maps of NGC 6951 (Kohno et al. 1999). Inside the ring 
(50 pc<r<350 pc), however, gas kinematics are compatible 
with regular rotation (Fig.|4j?). Of particular note, this is in clear 
contrast to the case of NGC 4579. The kinematics of the G- 
component, although compatible with circular motions, cannot 
be studied in detail due to insufficient spatial resolution. 

While most of the 3xlO 8 M molecular gas disk is feeding a 
starburst episode in the nuclear spiral arms at r~350 pc, a small 
amount of molecular gas (~a few 1O 6 M ) has been detected on 
the central engine. This central component reveals a prior ac- 
cretion episode down to scales of r~50 pc. It is unclear whether 
molecular gas detected in the bridge component is falling into 
the nucleus or migrating outward: the gas flow is fully compat- 
ible with regular circular motions from r=50 pc to r=350 pc. 

4. Gravitational torques and AGN fueling 

To explore more precisely the efficiency of feeding, we have 
estimated the gravitational torques exerted by the stellar poten- 
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tials (derived from the NIR images) on their molecular circum- 
nuclear disks (as given by the NUGA CO maps). After estimat- 
ing the role of stellar gravitational torques, we will investigate 
whether other mechanisms are required to explain the low level 
of nuclear activity in these galaxies. 

We first explain the general methodology employed and the 
basic assumptions in Sect. 14.11 The different steps are described 
in detail in Sect. 14. 1.1 1 and Sect. 14. 1.21 We discuss in Sect. IQ 
the results obtained from the application of this procedure to 
the four NUGA targets examined in this work. 

4.1. General methodology 

Gravitational forces are computed at each location in the plane 
of the galaxy, using near-infrared images to derive the underly- 
ing gravitational potential. We assume that the total mass bud- 
get is dominated by the stellar contribution and thus neglect 
the effect of gas self-gravity. We also assume a constant M/L 
ratio, and determine its best value by fitting the rotation curve 
constrained by the CO observations. From the 2D force field 
(Fx,Fy) we derive the torques per unit mass at each location 
(t(x,y)-x F y - y F x ). This torque field, by definition, is inde- 
pendent of the present gas distribution in the plane. The cru- 
cial step consists of using the torque field to derive the angular 
momentum variations and the associated flow time-scales. As 
explained below, the link is made through the observed distri- 
bution of the gas. 

With this aim, we assume that the measured gas column 
density (N(x, y)) derived from a CO intensity map at each offset 
in the galaxy plane is a fair estimate of the probability of find- 
ing gas at this location at present. In this statistical approach, 
we implicitly average over all possible orbits of gaseous par- 
ticles and take into account the time spent by the gas clouds 
along the orbit paths. We assume that CO is a good tracer of 
the total gas column density, since HI mass in the nuclei of 
galaxies is typically a very small fraction of the total gas mass. 
The torque field is then weighted by N(x,y) at each location 
to derive the time derivative of the local angular momentum 
surface density dL s {x,y) I dt-N(x,y) x t(x,y). 

In order to estimate the gas flows induced by these angu- 
lar momentum variations we produce azimuthal averages of 
dL s (x,y)/dt at each radius. The azimuthal average at each ra- 
dius, using N{x, y) as the actual weighting function, represents 
the global variation of the specific gas angular momentum oc- 
curring at this radius (dL/dt \g). Finally, the time-scales for gas 
inflow/outflow can be derived by estimating the average frac- 
tion of angular momentum transferred in one rotation. 

The validity of our estimate of the efficiency of stellar grav- 
ity torques to drive angular momentum transfer in the gas is 
based on the following simple hypothesis: we assume that the 
gas response to the stellar potential is roughly stationary with 
respect to the potential reference frame during a few rotation 
periods. We would like to stress that even in the particular case 
of nuclear bars, which might decouple from the outer stellar 
bars under certain circumstances, our assumption is still valid. 
When there are several stellar pattern speeds at different radii 
in a galaxy disk, numerical simulations show that the gas re- 



sponse tends to be coupled with the stellar potential pattern 
(Friedli & Benz [T9951 Garcfa-Burillo et al. [T9981 Bournaud 
& Combes 2002 2005 in prep.). The gas response adjusts its 
pattern speed to that of the dominant stellar pattern at a given 
radius. 

A different case is represented by non-axisymmetric per- 
turbations which can be driven by gas self-gravity and that are 
partly independent of or possibly decoupled from the stellar 
perturbations of the disk. Nuclear galaxy disks with a high gas 
surface density and a mostly axisymmetric stellar potential can 
be prone to develop this kind of gas self-gravitating pertur- 
bation. In this limiting case (not contemplated here), and al- 
though our calculation is still formally correct, the inclusion of 
gas self-gravity is required to derive the correct torque budget. 
In particular, if the gas disk decouples from the stellar pattern, 
the azimuthally averaged gravity torques exerted by the stellar 
pattern on the gas will very likely be close to zero. The main 
source for the torques if any should come from the gas instabil- 
ity itself. 

Note also that any radial variation of the CO-to-H2 con- 
version factor (Xco) is not expected to affect the estimated 
time-scales as these are independent of the global normaliza- 
tion factor as a function of radius. If Xco varies as a function 
of azimuth at a fixed radius, there could be a potential bias. 
However, observational evidence indicates that the dominant 
variation of Xco m the central regions of galaxies is radial 
(Solomon & Barrett fTWTl Regan et al-HSD- 

In the approach followed to derive the gravitational poten- 
tial we have assumed that the bulge is as flattened as the disk, 
and thus no attempt has been made to separate the bulge from 
the pure disk contribution. The implicit assumption of a highly 
flattened bulge is probably not wrong for some barred galaxies 
but for others it will overestimate the radial forces by at most 
a factor of 2 (e.g., Buta & Block l200fl l. On the other hand, the 
effect of bulge stretching due to deprojection can enhance the 
strength of bars, especially if they are aligned with the minor 
axis and the inclination angles are large. As this is not the case 
for the galaxies analyzed here, we instead expect that the value 
derived for the gravity torques will be typically underestimated 
in our case by a factor 1.5-2 although they will still have the 
same sign. 

4.1 .1 . Evaluation of stellar potentials 

The first step is to derive the stellar potential in the nuclear 
disks of these galaxies, using the high-resolution NIR images 
described in Sect. 12.21 Our working hypothesis is that NIR im- 
ages are less affected than optical ones by dust extinction or 
by stellar population biases (Quillen et al. 1994). The images 
are first deprojected according to the angles PA and i given 
in Tabled The images are then completed in the vertical di- 
mension by assuming an isothermal plane model with a con- 
stant scale height, equal to ~ 1/1 2th of the radial scale-length of 
the image. The potential is then derived by a Fourier transform 
method. We also assumed a constant mass-to-light (M/L) ra- 
tio, obtained by fitting the observed CO rotation curve-v, ,-for 
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each galaxy. The potential-<t>(/?, 6)- is then decomposed in the 
different m-modes: 



6) = ® (R) + V O m («) cos(m0 - cf> m (R)) 



(1) 



where 0> m (R) and <p,„(R) represent the amplitude and phase of 
the m-mode, respectively. 



Following Combes and Sanders d!981t . we define the 
strength of the ra-Fourier component, Q„,(R) as 



Q m (R) = mO m /R\F m 



(2) 



The corresponding strength of the total non-axisymmetric per- 
turbation is defined by: 



Qt(R) 



F max f ™ I ( BUM ) 

f T (K) r\ oe )„ 



F (R) 



dtbo(R) 
dR 



(3) 



where F'j' x (R) represents the maximum amplitude of the tan- 
gential force over all 8 and Fq(R) is the mean axisymmetric 
radial force. 

Figs. [S^fr to [HJa, Z? illustrate the quantitative description of 
gravitational potentials given by [Q,=i,2, Qr, <f>i=\,i] for the four 
galaxies examined in this work. 



angular momentum, respectively. More precisely, we evaluate 
the AGN feeding efficiency by deriving the average fraction of 
the gas specific angular momentum transferred in one rotation 
(Trot) by the stellar potential, as a function of radius, i.e., by the 
non-dimensional function AL/L defined as: 



AL 



dL 
dt 



1 

x - 
e L 



xr, 



t(R) 



xr r , 



(6) 



where Lg is assumed to be well represented by its axisymmetric 
average, i.e.,Lg = R x v rot . The absolute value of Lj AL deter- 
mines how long will it take for the stellar potential to transfer 
the equivalent of the total gas angular momentum. Assuming 
that the gas response to the stellar potential is stationary with 
respect to the potential reference frame during a few rotation 
periods, a small value of AL/L implies that the stellar potential 
is inefficient at present. The t(R) and AL/L curves derived from 
the 1-0 maps of the analyzed galaxies are displayed in Figs.|9j? 

toint>. 

To calculate how much gas mass is involved in the trans- 
fer driven by the stellar potential we have estimated the radial 
trend for the mass inflow (- sign)/outflow(+ sign) rate of gas 
per unit length as a function of radius (in units of M yr~'pc~' 
in Figs.ll3h to!16h) as follows: 



4.1 .2. Efficiency of gravitational torques 

After having calculated the forces per unit mass (F x and F y ) 
from the derivatives of <&(R, 6) at each pixel, the torques per 
unit mass-f(x, y)-can be computed by: 



t(x,y) = xF y -yF x 



(4) 



The sense of the circulation of the gas in the galaxy plane 
determines the sign of t(x,y): positive (negative) if the torque 
accelerates (decelerates) the gas at (x,y). We have next ob- 
tained the gravitational torque maps weighted by the gas col- 
umn densities derived from the CO 1-0 and 2-1 lines, N(x,y), 
i.e., we derive t(x,y) x N(x,y). These represent the effective 
variations of angular momentum density in the galaxy plane. 
We show in Fig. |9j? to 112b the normalized version of these 
maps, i.e., divided by [| N(x,y) x t(x,y) \] max for the 2-1 line 
(except for NGC 4321, where we used the 1-0 line). To esti- 
mate the radial gas flow induced by the torques, we have first 
computed the torque per unit mass averaged over the azimuth, 
using N(x,y) as the actual weighting function,i.e.: 



t(R) 



J e N(x,y)x(xF y -yF x ) 



(5) 



Results obtained for t(R) based on the CO(l-O) and CO(2- 
1) maps were seen to be virtually identical within the errors. 
Fig- EI 7 to!12t> show the results derived from the 1-0 line maps. 
By definition, t(R) represents the time derivative of the spe- 
cific angular momentum-L-of the gas averaged azimuthally, 
i.e., t{R)-dL/dt \g. To derive azimuthal averages, we assume a 
radial binning (AR) which corresponds to the original resolu- 
tion of the NIR images. Similarly to the torque maps, the sign 
of t{R), either + or -, defines whether the gas may gain or lose 



d 2 M 
dRdt 



dL 
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(7) 



where N(x, y) \ g is the radial profile of N(x,y) averaged over the 
azimuth for a radial binning AR. 

The inflow/outflow rates integrated out to a certain radius 
R can be derived as: 



dM d 2 M 



(8) 



dt dRdt 

Figs. 1 1 3Ej to 1 1 6t> display these integrated rates in units of 
Moyr" 1 . 

4.2. Tracking down gravitational torques in NUGA 
targets 

4.2.1. NGC 4321 

Fig- 13* shows the K-band image of the nucleus of NGC 4321, 



obtained by Knapen et al. ( 1995 1, deprojected onto the galaxy 
plane (Xc/Yc coordinates). The NIR image shows a nuclear 
bar (denoted as BAR(n) in Fig.|5Jj) which is roughly parallel to 
the large-scale 8 kpc stellar bar (oriented along PA=1 12° in sky 
coordinates). The nuclear bar strength is the main contributor 
to Q r , which is typically -0.18 for r<500 pc (Fig.|5j?). 

Fig. |9jj, shows that the derived torques change sign fol- 
lowing a characteristic 2D butterfly pattern. Quadrants I-to-IV 
will define hereafter the regions where the signs of the torques 
driven by the dominant perturbation of the stellar potential are 
expected to be constant (Figs l9H2"h ). Assuming that the gas 
circulation is counterclockwise, Fig. |9jj shows that the bulk 
of the CO emission along the spiral arms lies at the trailing 
edges of the nuclear bar where torques are positive (quadrants 
I(+) and !!!(+)). While there is significant CO emission arising 
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NGC4321 
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Fig. 9. a)(left) We overlay the 12 CO(1-0) contours with the map of the effective angular momentum variation (t(x,y) x N(x,y), 
as defined in the text) in the nucleus of NGC 4321. The grey scale is normalized to the maximum absolute value in the map, 
i.e., [| N(x,y) x t(x,y) \] max . The derived torques change sign as expected if the butterfly diagram, defined by the orientation of 
quadrants I-to-IV, can be attributed to the action of the nuclear bar of NGC 432 1 . b)(right) The torque per unit mass averaged over 
azimuth-f(r)-and the fraction of the angular momentum transferred from/to the gas in one rotation-c/L/L-are plotted. Torques 
are strong and positive for the bulk of the molecular gas in NGC 432 1 . This includes the vicinity of the AGN. Torques are negative 
but comparatively weaker on intermediate scales (r=250-550 pc) and in the outer disk r>900 pc. 



from the leading quadrants of the bar (quadrants II(-) and IV(- 
)) both on intermediate scales (r=250-550 pc) and in the outer 
disk r>900 pc, negative torques there are weaker than posi- 
tive torques measured over the trailing quadrants (see Fig.[5f>). 
Most remarkably, closer to the AGN (r<200 pc), the sign of the 
average torques t{R) is also positive. Although the central CO 
source is marginally resolved in the CO map of Garcfa-Burillo 
et al. ( 1998 ), and, therefore, higher resolution observations are 
required to settle the question, positive torques seem to domi- 
nate over negative ones at r<200 pc. 

In terms of mass inflow rates, the overall feeding budget in 
NGC 4321 is clearly positive at all radii as shown in Fig. 1131 
This would imply that the radial gas flow goes outward un- 
less other mechanisms more efficiently transfer the gas angu- 
lar momentum outwards. At r~100 pc, dM/dt~+0.4Mo yr _1 
and the radial rate is still positive at r~500 pc, where 
dM/dt~+O.4M0 yr" 1 . Similarly to the case of the other tran- 
sition object, NGC 4826 (see Sect. WF2b . AGN feeding 
driven by stellar torques is quenched at present in NGC 4321. 
Although the gas concentration is large close to the AGN, the 
present gas configuration favors feeding only on intermediate 
scales (r=250-550 pc). This result is compatible with the CO- 
based diagnostic of the gas kinematics discussed in Sect l3.ll 



4.2.2. NGC 4826 

Fig-Eh shows a zoom into the central ~200 pc region of the de- 
projected H-band HST image of NGC 4826. We plot in Fig.gJ? 
Q,( z'=l,2), (pi (i=l,2) and inside the full field-of-view of 
the image (r=530 pc). As shown in Fig.[|J>, the deviations from 
axisymmetry of the stellar potential are exceedingly small at all 
radii (Qj ~0. 02-0. 05) except for r<75 pc where a weak oval 
perturbation is detected. The latter is denoted as OVAL(n) in 
Fig.|6jj. This oval seems slightly off-center with respect to the 
AGN; this would account for the contribution from am- 1 
stellar mode to Q^, identified in Fig.[6j>. 

As already suggested by the first calculations of paper I, 
we find that the weak stellar perturbations present in the in- 
ner 500 pc of NGC 4826 are inefficient at feeding the AGN 
at present. Fig. HOfe shows how the derived torques change 
sign in the central 200 pc of the disk. The distribution of 
torques, defined by quadrants I-to-IV in Fig.llOb. shows the ex- 
pected behavior if torques are mostly due to the oval distortion 
of NGC 4826. Moreover, stellar torques are weak and mostly 
positive inside r~150 pc (Fig. \10b). Very close to the AGN 
(r<20 pc), the sign of average torques might become marginally 
negative, but as we lack spatial resolution on these scales, the 
torque estimate becomes uncertain. 
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Fig. 10. Same as Fig.[9]but for NGC 4826, using the I2 CO(2-l) map in the overlay and the I2 CO(1-0) map in the radial averages. 
The derived torques change sign as expected if the butterfly diagram, defined by the orientation of quadrants I-to-IV, can be 
attributed to the action of the oval distortion of NGC 4826. Stellar torques are weak and mostly positive in the nucleus of 
NGC 4826 at present. 

NGC4579 
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Fig. 11. Same as Fig.llOlbut for NGC 4579. The torques change sign as expected if the butterfly diagram, defined by the orientation 
of quadrants I-to-IV, can be attributed to the action of the large-scale bar of NGC 4579. Torques are systematically strong and 
negative for the bulk of the molecular gas in NGC 4579, from r=200 pc out to r=1200 pc. In the vicinity of the AGN, however, 
torques become positive and AGN feeding is not presently favored. 
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Fig. 12. Same as Fig.llOlbut for NGC 695 1 . The torques change sign as expected if the butterfly diagram, defined by the orientation 
of quadrants I-to-IV, can be attributed to the action of the nuclear oval rather than to that of the outer bar of NGC 695 1 . Torques 
are systematically negative over the nuclear spiral arms down to the inner radius of the pseudo-ring, r~300 pc. In the vicinity of 
the AGN, torques become positive and AGN feeding is not presently favored. 



Alternatively, the lopsided gas disk may be unrelated to the 
stellar oval. Instead, gas self-gravity could be driving the gas 
response in the nucleus NGC 4826 (see discussion in paper 
I). In this case the gas disk may be decoupled from the stellar 
pattern, and therefore the latter would be viewed as a mostly 
axisymmetric component from the system of reference of the 
gas, i.e., providing zero torque on the gas disk. The main source 
for the torques should come from the gas instability itself. In 
NGC 4826 m = 1 instabilities do not seem to favour gas inflow, 
however (paper I). In either case, this is compatible with the 
scenario depicted in Sect. 13.21 in which there is currently no 
evidence of fueling of the AGN. 

In spite of the large molecular gas concentration close to 
the AGN (r<150 pc), the overall feeding budget in NGC 4826, 
quantified by the mass inflow rates shown in Fig. [21 is clearly 
positive inside the full field-of-view of the NIR image. This 
implies that the predicted radial mass flow should be dom- 
inated by outflow rather than inflow motions. At r~100 pc, 
dM/dt~+O.2M0 yr~' and the inflow rate is still positive at 
r~500 pc, where dM/dt~+O.45M yr 1 . 

4.2.3. NGC 4579 

Fig-Et shows the central r~l kpc region of the I-band HST 
image of NGC 4579 deprojected onto the plane of the galaxy. 
The dominant stellar perturbation in the HST image of the disk 



is the large-scale 9 kpc bar of NGC 4579, denoted BAR and 
oriented as shown in Fig.0;. As illustrated by Fig.^t?, the po- 
tential strength of the bar is dominated by an m - 1 mode for 
r<250 pc where Q r >0.10. 

Fig.lllfe shows the 2D pattern of gravitational torques for 
r<l kpc. Stellar torques change sign as expected if the orien- 
tation of quadrants I-to-IV can be attributed to the large-scale 
stellar bar of NGC 4579. Assuming that the sense of gas cir- 
culation is clockwise in the galactic plane, Fig. II lb shows that 
the bulk of the CO emission along the spiral arms lies where 
torques are strong and negative. Fig II lb shows that the az- 
imuthally averaged torques are strong and negative from the 
edge of the image at r=1200 pc (and most probably from larger 
radii) down to r=200 pc. 

Closer to the AGN (r<200 pc), the bulk of the CO emis- 
sion is concentrated North of the nucleus in the trailing quad- 
rant I(+), where stellar torques are positive. The CO complex 
closest to the AGN is in the leading quadrant II(-), East of 
the nucleus at r~150pc, where torques are negative. However, 
Fig. II lb shows that the azimuthally averaged torques are pos- 
itive for r<200 pc. As discussed in Sect. 13. 31 the interpretation 
of the complex kinematics observed for molecular gas in this 
region is not straightforward but does not favor a direct feeding 
of the AGN. One of the possible scenarios explored by Garcfa- 
Burillo et al. 2005 (in prep.) invokes that gas is being entrained 
by expanding motions at r<200 pc. Note that in this case we do 
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not expect the gas pattern to be stationary with respect to the 
stellar pattern and therefore, the torque calculation developed 
above may not be relevant on scales r<200 pc. Nonetheless 
this would not change the picture where the central engine of 
NGC 4579 is not being fueled at present. 

In summary, the overall mass inflow budget is clearly neg- 
ative down to r~200 pc due to the action of the large-scale bar 
of NGC 4579 (Fig. [BJ. The radial flow expected for the gas 
is clearly inward down to r~200 pc. At r~500 pc, dM/dt — 

0. 6Mq yr . Inside r~200 pc stellar torques do not currently 
help feeding. 

4.2.4. NGC 6951 

Fig. [SJ; shows the deprojected J-band HST image of the nu- 
clear region (r<900 pc) of NGC 6951. The HST image allows 
the identification of the central 1.8 kpc of the large-scale bar 
of NGC 6951 (denoted as BAR in Fig.EJj). The bar clearly ex- 
tends beyond the image field-of-view (r=900 pc). Moreover we 
report on the detection of a small oval distortion most promi- 
nent inside r=400 pc (denoted as OVAL(n) in Fig. |8ji). The 
m — 2 mode of the oval perturbation dominates the strength of 
the potential at r<400 pc. While it is true that we likely under- 
estimate the contribution of the large-scale bar to the strength 
of the m — 2 mode at r<400 pc due to the restricted field-of- 
view of the HST image, we do not expect this contribution to be 
significant on these scales, considering the low Ch values typ- 
ically shown by large-scale bars at radii of a few 100 pc. The 
measured strength of the oval is moderate-to-low: Qj- ~0.06, 

1. e., quite comparable to the strength of the oval potential de- 
tected in the nucleus of NGC 4826. We note that the strength 
of the large-scale bar is likely much higher than that of the oval 
but it is relevant only for the gas which lies beyond the r~900 pc 
circumnuclear CO disk. 

The NGC 6951 torque map shown in Fig. 112b suggests in- 
deed that it is the nuclear oval which shapes the NGC 6951 
torques map. The orientation of quadrants I-to-IV associated 
with the nuclear oval fits the pattern of the torques much better 
than that associated to the bar. As the sense of gas circulation 
is clockwise in NGC 6951, Fig. 112b shows that CO emission 
over the spiral arms in quadrants I(-) and III(-) lies at the lead- 
ing edges of the nuclear oval where torques are negative down 
to the inner radius of the pseudo-ring, at r~300 pc. While the 
torques due to the large-scale bar on the inner CO ring are weak 
at present, torques are expected to be strong and negative at 
the leading edges of the bar. The offset dust lanes detected in 
the optical images of NGC 6951 (e.g., Perez et al. 120001 in- 
dicate that gas inflow is at work along the bar on large scales. 
The combination of the bar and the nuclear oval has created a 
pattern of torques which are negative down to r~300 pc. The 
bulk of the molecular gas in the central 5 kpc of NGC 695 1 
has already fallen from the xi orbits down to the X2 orbits of 
the bar where it feeds a nuclear starburst in the ring. In the 
vicinity of the AGN, however, the average stellar torques be- 
come marginally positive and AGN feeding is not favored by 
the present configuration of the stellar potential. 



Fig. ^] shows the overall mass inflow budget for 
NGC 695 1 . The radial flow expected for the gas is clearly in- 
ward down to r~300 pc. At r~500 pc, dM/dt— 0.25M Q yr -1 . 
Inside r~300 pc, however, stellar torques do not favour AGN 
feeding at present in this Seyfert. 

5. Discussion and conclusions 

The first results obtained from the analysis of stellar torques 
in the circumnuclear disks of a subset of NUGA targets re- 
veal complex trends in the details of the feeding history along 
the Transition/LINER/Seyfert sequence. The overall picture 
emerging in the two Seyferts/LINERs (i.e., NGC 6951 and 
NGC 4579) suggests that AGN feeding may proceed in two 
steps. In a first step, gravity torques do help, by getting rid of 
gas angular momentum, driving the gas inwards and feeding a 
nuclear starburst on scales of ~a few 100 pc. On smaller scales, 
however, stellar torques play no role in AGN fueling in the cur- 
rent epoch: torques on the gas are not negative all the way to the 
center, but on the contrary they become positive and quench the 
feeding. This may explain why molecular gas seems to 'avoid' 
the inner 200-300 pc of NGC 6951 and NGC 4579 where we 
measured M gas <a fewxl0 5 -10 6 MQ (see also the properties of 
NGC 7217 as discussed by Combes et al. 2004i. Observational 
evidence indicates that in NGC 4579 molecular gas on these 
scales is presently flowing outward. In a second step, a mech- 
anism complementary to gravity torques may be required to 
drive gas inflow on smaller scales. 

Most notably, a much larger molecular gas concentration 
has been driven inwards in the transition objects analyzed in 
this paper: ~3xlO 7 -lO 8 M0 inside r=150 pc. However, we 
have found that stellar torques are also unable to drain angu- 
lar momentum from the massive circumnuclear gas disks of 
NGC 4321 and NGC 4826. The underlying mechanism respon- 
sible for halting gas inflow inside r=150 pc seems to be differ- 
ent in these two transition objects. On similar spatial scales, 
the stellar torques created by the nuclear bar of NGC 4321 are 
much stronger than those created by the weak oval perturbation 
of NGC 4826. Furthermore, the m = 1 perturbation identified 
in the gas disk of NGC 4826, which is probably unrelated to 
the oval perturbation, could be the main donor of angular mo- 
mentum to the gas at these radii (paper I). 

Our quest for a mechanism of AGN feeding based only 
on stellar torques has led to a paradoxical result for the four 
cases studied here: inflow is actually thwarted close to the 
AGNs (r<100-200 pc). One possible explanation is that the re- 
sponsible agent in the stellar potential could be transient, or 
as short lived as an individual AGN episode (hence: <10 7 8 
years). The feeding phase from 100 pc to 1 pc could be so 
short that the smoking gun evidence in the potential is likely 
to be missed. Alternatively, this temporary inability of gravi- 
tational torques could be overcome by other mechanisms that, 
over time, become competitive with non-axisymmetric pertur- 
bations. Among the different mechanisms usually cited in the 
literature, dynamical friction and viscous torques have been in- 
voked to help AGN feeding on ~100 pc scales. Gravitational 
torques can become positive close to the central engine as il- 
lustrated above, in contrast with viscous torques or dynamical 
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friction that always favour gas inflow. This implies that stellar 
torques could regulate the gas flows in galactic nuclei in com- 
bination with other mechanisms that secularly drain the gas an- 
gular momentum. 

We discuss below the efficiency of dynamical friction and 
viscosity versus gravity torques for driving AGN feeding in 
the four galaxies analyzed in this paper. This efficiency can be 
quantitatively estimated by comparing the typical time-scales 
of these processes with those of the gravitational torques de- 
rived above for each galaxy. 



Galaxy 




T 


Tdf 


T 

1 vise 






(Trot) 


(Trot) 


(Jrot) 


NGC 4321(r= 


=200 pc) 


8 


30 


20 


NGC 4826(r= 


=50 pc) 


70 


5 


5 


NGC 4579(r= 


=200 pc) 


8 


30 


22 


NGC 6951(r= 


=200 pc) 


50-100 


30 


20 



Table 2. We compare the characteristic time-scales of gravity 
torques (T grav ), dynamical friction (Tdf) and viscosity (T,, iJC ) in 
the nuclear disks of NGC 4321 , NGC 4826, NGC 4579 and 
NGC 695 1 . Time-scales are given in units of rotation periods 

{T ro t). 



5.1. Assisting gravity torques 
5.1.1. Dynamical friction 

Dynamical friction of giant molecular clouds (GMCs) in 
the stellar bulge of a galaxy is often invoked as a possible 
mechanism of fueling AGN. According to the well-known 
Chandrasekhar (1943) formula, the time-scale of dynamical 
friction for a GMC of mass Mqmc at a radius r where rota- 
tional velocity is V m would be T df ~10 8 yr (r/100 pc) 2 x(10 6 
M /M GMC ) _1 x(y ro ,/2OOkms _1 ). We have derived the decay 
time-scale (T d f) of a typical GMC (M CM c ~1O 6 M ) in the nu- 
clear disks of the 4 galaxies discussed above. We derive Tdf 
at the radii where we have identified a barrier of positive grav- 
ity torques and get V rot from the CO observations. The values 
of T,/y are typically ~30, except in NGC 4826 where we ob- 
tain formally Tdf ~5 (see Table |2jl. The comparison between 
T d f and the corresponding time-scales for gravity torques T gmv 
(defined as T gmv ~ L/AL from Eq. |6)) implies that friction is a 
relatively slow process compared to gravity torques, except in 
NGC 4826 and NGC 6951. 

Despite these results, we suspect that the values estimated 
for Tdf are a strict lower limit and so severely overestimate the 
efficiency of dynamical friction in all cases. The reason for this 
assessment is that Chandrasekhar's formula implicitly assumes 
that a GMC is a bound rigid body, with a highly concentrated 
mass. This approximation underestimates the real value of Tdf. 
Molecular clouds are indeed a loose ensemble of clumps of 
different sizes, resembling fractals. The typical masses of the 
densest fragments can be as low as 1CT 3 M G (e.g., Pfenniger 
& Combes [l994i . Therefore the scattering efficiency of a re- 
alistic GMC unit of 10-30 pc size including a complex mass 
spectrum of clumps inside is much less than estimated above: 
the Coulomb parameter of a GMC in Chandrasekhar's formula 
should tend to vanish. We therefore conclude that dynamical 
friction is likely to be quite a slow process, which to first ap- 
proximation can be neglected relative to gravity torques. 



5.1.2. Viscosity 

Viscous torques are generally weak, and their associated time- 
scales (T V i s ) are quite long at large radii. The reason behind 
this is that T,,, s grows as the square of the radius according to 
the classical diffusion equation of Pringle ( 1981 1 which gives 



the time evolution of an axisymmetric disk distribution, 2(r,f), 
under the action of viscous transport: 



at 



'1 

dr 



d[vT.r\dQ./dr)] 



dr 



'cKClr 2 )) 
dr 



(9) 



In this formula, the kinematic viscosity, v, is usually de- 
scribed as A x <t v , where A is the characteristic mean-free path 
for viscous transport and <x,, is the gas velocity dispersion. 
Eq. (9) allows to derive by simple dimensional analysis the typ- 
ical radial dependence of the time-scale for viscous transport in 
an axisymmetric disk, if we assume that the disk rotation curve 
can be described by a power law (i.e., V rof oc r l ~"): 



2 — a r 

'vis ~ X 

a v 



(10) 



However, as illustrated below, the exact value of T v , s 
strongly depends on the assumed initial density distribution law 
for the disk. Fig. 117b shows the time evolution of a smooth disk 
distribution under the action of viscous torques. In this calcu- 
lation we have taken y=100 pc km s _1 , assuming A-10 pc and 
cr,,=10 km s . After 200 Myr the radial gas flows are very 
small. However, the efficiency of viscous transport can be en- 
hanced if the initial distribution of the disk is characterized by 
strong density gradients. This is the case of nuclear contrasted 
rings, especially when these are located in the inner regions of 
galaxies (r~100-500pc) as frequently is the case for early-type 
barred spirals. Moreover, and contrary to the conventional wis- 
dom (based on low resolution observations) that most galaxy 
rotation curves are close to rigid body for r<500 pc, Q. can show 
a strong variation with radius in the inner regime. Therefore 
galactic shear can still be very high on these scales; this also 
favours viscous transport. 

Fig. I17t > illustrates the time evolution of a nuclear ring of 
r~100 pc under the action of viscous torques in a galaxy disk. 
The radial gas flows are remarkably faster than in the pre- 
vious case: in typically 5xl0 7 yr, the ring starts to dissolve 
and a significant amount of gas can reach the center of the 
galaxy. We compare in Table |2 the characteristic time-scales 
of viscous torques and gravity torques in the 4 galaxies an- 
alyzed in this paper. Power-law rotation curves have been fit 
within the relevant radial ranges adapted to each galaxy, as- 
suming in all cases v=100 pc km s , i.e., a viscosity pre- 
scription similar to the the axisymmetric smooth disk solution. 
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Viscosity seems a viable mechanism to drive significant gas 
inflow in NGC4826 at 1-50-100 pc, while it is probably in- 
efficient against the stronger gravity torques in NGC 4579 and 
NGC 432 1 for r<200 pc. Moreover, gas lying inside the nuclear 
ring of NGC 695 1 could be falling to the AGN favoured by vis- 
cous torques at r<200 pc. Note that these results are unchanged 
if we allow gravity torques to be a factor of 1 .5-2 stronger, as 
argued in Sect. 14. II 

We can conclude that, in spite of all the uncertainties in the 
derivation of time-scales, viscous torques are the most viable 
mechanism to generate gas inflow on scales ~ 100-200 pc if 
they act on a contrasted nuclear ring distribution, and do not 
have to fight against very strong positive torques from gravity. 

5.2. A scenario for self-regulated activity in LLAGNs 

In this work we have shown that gravity torques exerted by the 
stellar potential on the gas disks of four LLAGNs, purposely 
chosen to represent the whole range of activity classes in the 
NUGA sample, fail to account for the feeding of the AGN at 
present. We conclude that gravity torques need to be assisted 
in due time to drive the gas to the center. Based on these re- 
sults, we develop in this section a simplified general scenario 
in which the onset of nuclear activity can be understood as a re- 
current phase during the typical lifetime of any galaxy. In this 
scenario the recurrence of activity in galaxies is indirectly re- 
lated to that of the bar instabilities although, as argued below, 
the active phases are not necessarily coincident with the max- 
imum strength of a single bar episode. While bars can build 
up gas reservoirs towards the central regions of galaxies in 
the shape of nuclear rings in early type objects, the dynami- 
cal feed-back associated with these gas flows can destroy or 
considerably weaken the bars (Norman et al. 119961 Bournaud 
& Combes 2002 1. When this happens, gravity torques are neg- 
ligible and can make way for other competing mechanisms of 
gas transport, such as viscous torques, in an almost axisymmet- 
ric system. 

The first step in this evolutionary scenario could begin with 
an axisymmetric disk of gas. Shortly after the galaxy disk 
would be prone to a bar instability. If the galaxy is early-type, 
i.e., with a sufficient bulge or central mass concentration, then 
the bar pattern is likely to have one or even two Inner Lindblad 
Resonances (ILRs). A nuclear ring would form either close to 
the single ILR or between the two ILRs. At this stage, the gas 
would be driven inwards from corotation to the ILR, and would 
accumulate in the ring where a burst of star formation can oc- 
cur. However positive gravity torques would prevent the gas 
from flowing further in: under the bar forcing, gas inside the 
ring could be evacuated outwards (Combes 1988 1. The bar thus 
would restrict the amount of fuel available for AGN feeding to 
that lying very close to the center where it is under the domi- 
nant gravitational influence of the black hole. As the potential 
would be mostly axisymmetric in the central ~10 pc, gas could 
remain there feeding the central engine. 

The infall of gas driven by a bar is self-destructive, how- 
ever: it progressively weakens and destroys the bar, so that 
the potential returns to axisymmetry (Norman et al. 119961 



Bournaud & Combes 2002 1. At this stage gas piled up in the 
nuclear ring could dissolve to form a smoother disk through 
viscosity which may have become competitive against gravity 
torques. The infall of gas produced by viscous torques could 
replenish the exhausted fuel near the central engine triggering 
a new phase of nuclear activity. 

In this scenario, the feeding of active nuclei in early-type 
spirals is contemplated as a two step process: first gravity 
torques bring the gas of the large-scale disk to the nuclear 
ring, and when the dynamical feedback has destroyed the bar, 
the viscous torques could smooth out the ring, and bring gas 
into the central 10 pc, where it is under the influence of the 
Keplerian potential of the black hole. The disk becomes ax- 
isymmetric and the cycle can be restarted at its first step. In 
particular, the disk will be prone to a new bar instability if gas 
is accreted from the outer parts of the disk. These multiple bar 
phases must occur in gas-rich spiral galaxies, in order to ac- 
count for the observed frequency of bars (Block et al. 2002; 
Bournaud & Combes l2002ll2004b . 

We relate the recurrence of episodes of bar formation and 
destruction with the fueling of active nuclei in galaxies. The 
self-regulated competition between gravity torques and viscos- 
ity in galaxy nuclei may lead typically to several episodes of 
activity (each lasting for ~ 10 8 yr) during a single cycle of 
bar formation/destruction (lasting for 10 9 yr). These activity 
episodes are not expected to be strongly correlated with the 
phases of maximum strength for the bar, but they may appear 
at different evolutionary stages of the bar potential, depending 
on the balance between gravity torques and viscosity. The ob- 
served prevalence of outer rings in Seyferts, interpreted as a 
sign of bar dissolution, would support this scheme (Hunt and 
Malkan 1999 1. Furthermore, the ample variety of morphologies 
revealed by the CO maps of the circumnuclear disks of NUGA 
targets corroborates that there is no universal pattern associated 
with LLAGNs (Garcfa-Burillo et al. 2004 1. Activity can also be 
found in a galaxy during its axisymmetric phase (e.g., see the 
case of the NUGA galaxy NGC 5953, discussed by Combes et 
al. 2005, in prep). 

The gas directly responsible for the Seyfert activity in 
galaxies like NGC 4579 or NGC 695 1 must have been brought 
to the center during the previous axisymmetric phase, while the 
bar at present is emptying the region inside the nuclear ring, 
thus regulating the amount of gas available for the active nu- 
cleus. In the case of NGC 4579, gravity torques clearly over- 
come viscosity, while in NGC 695 1 the balance between both 
mechanisms might be reached soon. In the case of NGC 4826, 
the nucleus has a very low level activity in spite of the presence 
of a large amount of gas in the vicinity, since the potential is 
almost axisymmetric, and the nuclear ring could be undergo- 
ing dissolution by viscosity. The active phase would develop 
during the next phase, while the disk would soon become bar 
unstable again. However, the role of the m = 1 gas perturba- 
tions identified in the disk of NGC 4826 might be to slow down 
the dissolution of the ring, as kinematics still suggest outward 
motions close to the AGN. The case of NGC 4321 is still early 
in the time evolution, as the bar is only now entering the de- 
struction phase, through the formation of a nuclear bar inside 
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the resonant ring. This nuclear bar efficiently prevents nuclear 
feeding through gravity torques. 

The observed fraction of LLAGNs in the Local Universe 
(-44%, including LINERs; Ho et al. H997i could impose tight 
constraints on the expected number of activity episodes per bar 
formation/destruction cycle in this scenario. Based on these 
statistics, a rough estimate would raise this number to ~4. 
However this estimate can only be taken as mostly specula- 
tive for the time being. The principal source of uncertainty, 
which prevents us from making a more quantitative predic- 
tion, resides in the fact that there is no consensus on what is 
the typical duration of a nominal AGN duty-cycle. The differ- 
ent values of AGN duty-cycles discussed in the literature range 
from 10 7 to 10 8 yr (Ho et al. 120051 Martini 120031 Wada|2004 
Merloni 2004). Moreover, the required number of individual 
episodes could be lowered to ~a few if we exclude starburst 
dominated LINERs from the LLAGN population or even more 
if we restrict the activity classification to Seyferts, the most ac- 
tive members in the family. 

In barred galaxies with comparatively less prominent 
bulges, i.e., characterized by the absence of an ILR barrier, vis- 
cosity may not be needed to take over gravity torques close to 
the AGN. The gas flow should not be stalled at a radius of a 
few ~100 pc forming a nuclear ring in this case. There are a 
few examples in the NUGA survey of strongly barred galax- 
ies showing a strong nuclear concentration of gas instead of a 
ring. This class of ILR-free bars could be more common in the 
early Universe during the formation of disks and the growth of 
supermassive black holes (Sellwood & Shen 2004 1. 

Furthermore, other mechanisms different from axisymmet- 
ric viscous flow could be efficient at overcoming the ILR bar- 
rier imposed by most stellar bars in the Local Universe. As 
mentioned above, the ultimate agent for the feeding could still 
reside in the stellar potential but would be transient and their 
detection might be elusive. In addition, the role of nuclear 
m — 1 and m = 2 gaseous spirals in AGN feeding is still not 
completely elucidated. Most interestingly, a few examples of 
these gas instabilities are found in the available CO maps of 
AGNs, including some of the galaxies of the NUGA sample 
(Garcfa-Burillo et al. 120031 . 

Finally, the role of gas self-gravity to drive gas inflow 
may not be negligible in those cases where the distribution 
of gas is a significant source of non-axisymmetry for the total 
gravitational potential, especially if the stellar potential itself 
appears as featureless or mostly axisymmetric. We can esti- 
mate the influence of gas self-gravity by calculating the typi- 
cal gas mass fractions as a function of radius in the case of the 
galaxies studied in this paper. In NGC 6951 and NGC 4579, 
the gas mass fractions are very low throughout the disk from 
r~100 pc (-1%) up to r~1000 pc (-3%). Therefore we do not 
expect gas instabilities to have a significant influence in the 
total gravitational potential in these two galaxies. In contrast, 
gas mass fractions are comparatively larger at the same radii in 
NGC 4826 and NGC 432 1 : ~ 1 0- 1 5 % at r~ 1 00 pc in both galax- 
ies. This gas mass fraction is still high at r~1000 pc (>10%) 
in NGC 4321, while it decreases down to -5% in NGC 4826 
at that distance. This evidence points to a larger influence of 



gas self-gravity in NGC 4826 (as discussed in paper I) and 
NGC 4321 (as discussed by Wada et al. fT998l . 

These results highlight the need for a significant increase 
in the number objects for which high-resolution CO maps are 
available together with careful case-by-case studies in the quest 
for evidence of AGN feeding. The scenario proposed in this 
work, as well as any alternative model accounting for nuclear 
activity in galaxies, remains to be tested using a larger sample 
of LLAGNs. 
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Fig. 5. a) (upper panel) We overlaid the 12 CO(1-0) PdBI in- 
tensity map (contours) on the K-band image of Knapen et 
al. {1995 1 (grey scale) obtained for the nucleus of NGC 4321; 
both images have been deprojected onto the galaxy plane. Units 
on X/Y axes (AXc/AYc) correspond to arcsec offsets along the 
major/minor axes with respect to the AGN. The nuclear bar- 
BAR(n)-is parallel to the large-scale bar-BAR. b) (lower pan- 
els) Strengths (Q,, z'=l,2) and phases (</>,, z=l,2) of the m = 1 
and m = 2-Fourier components of the stellar potential inside 
the image field-of-view (r=1100 pc). The 0,-angles are mea- 
sured from the +X axis in the counter-clockwise direction. The 
total strength of the potential is represented by Qj. Q-values 
for r<50 pc (not shown) are blanked due to insufficient grid 
sampling close to the nucleus. 
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Fig. 6. a) (upper panel) We overlaid the 12 CO(2-l) PdBI in- 
tensity map (contours) on the H-band HST image (grey scale) 
obtained for the nucleus (r<200 pc) of NGC 4826; both images 
have been deprojected onto the galaxy plane. The orientation 
of the nuclear oval perturbation-OVAL(n)-is shown, b) (lower 
panels) As in Fig. 13 we plot (Q„ z'=l,2), (0,, z'=l,2) and Q T 
inside the image field-of-view (r=530 pc) with an inner trunca- 
tion radius r=10 pc. 
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Fig. 7. a) (upper panel) We overlaid the 12 CO(2-l) PdBI in- 
tensity map (contours) on the I-band HST image (grey scale) 
obtained for the nucleus of NGC 4579; both images have been 
deprojected onto the galaxy plane. The orientation of the large- 
scale 9 kpc bar-BAR-is shown, b) (lower panels) As in Fig. [5] 
we plot (Q,-, i=l, 2), ((/>[, i-1,2) and Q T inside r=1200 pc with 
an inner truncation radius r=50 pc. 
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Fig. 8. a) (upper panel) We overlaid the 12 CO(2-l) PdBI in- 
tensity map (contours) on the J-band HST image (grey scale) 
obtained for the nucleus of NGC 695 1 ; both images have been 
deprojected onto the galaxy plane. We indicate the orientation 
of the large-scale bar-BAR-and of the nuclear oval distortion- 
OVAL(n)-that is detected in the HST image at r<500 pc. b) 
(lower panels) As in Fig. we plot (Q,, z'=l,2), (<pj, z'=l,2) 
and Qt inside the image field-of-view (r=900 pc) with an inner 
truncation radius r=25 pc. 
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Fig. 13. a)(upper panel) We represent the radial variation of 
the mass infiow(-) or outfiow(+) rate of gas per unit radial 
length in the nucleus of NGC4321 due to the action of stellar 
gravitational torques. Units are Mq yr _1 pc _1 . b)(lower panel) 
Here we plot the mass inflow/outflow rate integrated inside a 
certain radius r in Mq yr _1 . As is the case of the other tran- 
sition object, NGC 4826, the overall budget in NGC4321 is 
clearly positive at all radii. 
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Fig. 14. Same as Fig.fllfc. b but for NGC 4826. As is the case 
of the other transition object, NGC 4321, the overall budget in 
NGC 4826 is clearly positive at all radii. 
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Fig. 15. Same as Fig. [jjj?, b but for NGC 4579. The overall 
mass inflow budget is clearly negative down to r=300 pc due 
to the action of the large-scale bar. Inside this radius, stellar 
torques do not favour AGN feeding in this LINER/Seyfert. 
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Fig. 16. Same as Fig. [l3j?, b but for NGC 6951. The overall 
mass inflow budget is clearly negative down to r=300 pc due to 
the combined action of the large-scale bar and the nuclear oval. 
Inside this radius, stellar torques do not favour AGN feeding in 
this Seyfert. 
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Fig. 17. a) The time evolution of gas surface density in an ax- 
isymmetric disk, where the rotation curve is of the form: V rof oc 
r , with a = 0.1. To normalize, we fit V ro ,=200 km s~ l at 
r-\ kpc. Initially the gas is distributed according to a surface 
density law proportional to £ oc 1 + cos(nr/r c ), until the cut-off 
radius r c = 0.3 kpc. We represent Logio£(r,t) in arbitrary units. 
The various curves are the evolution for times spaced by 2xl0 7 
yr, until 2xl0 8 . The column density profile shows hardly any 
change after 2xl0 8 yr. b) Same as a) but here for an initial gas 
distribution in the shape a narrow ring between 90 and 1 10 nc 



